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We review, from an experimental point of view, the current status of ultra-relativistic 

■ nuclear collisions with heavy beams. 

V ; 1. Introduction 

Reactions between heavy nuclei at ultra-relativistic energies have now been studied for 
a number of years at the BNL AGS (11.4 GeV/nucleon) and at the CERN SPS (158 

■ GeV/nucleon). Progress has been rapid but it is only in the last year that sufficiently 
complete and comprehensive data sets have become available that a "picture" is beginning 
to emerge. In the following we will briefly review the experimental highlights since the 
last Quark Matter conference at Heidelberg. Particular emphasis will be placed on the 
identification of possible collective behavior and hydrodynamic flow as well as on the 
question whether or not there is local equilibrium at freeze-out. As most progress since 
QM'96 in the leptonic sector will only be reported at this conference we will only touch 
upon recent experimental developments in the study of the low-mass dilepton continuum 
and very briefly summarize the status of anomalous charmonium suppression. 



2. Global Variables 

2.1. E T and dE T /dr] Distributions 

Measurements of the transverse energy production and its spatial distribution in cen- 
tral Au+Au collisions at the AGS |I[ and Pb+Pb collisions at the SPS |2j have yielded 
maximum pseudorapidity densities of dE^/dr? = 200 and 450 GeV, respectively. Using a 
simple Bjorken-type estimate these values imply energy densities in the fireball formed in 
the collision of about 1.3 and 3 GeV/fm 3 . Recent results from solutions of QCD on the 
lattice || imply critical temperatures (for systems including dynamical quarks) well below 
200 MeV. The corresponding critical energy density is then of the order of 1-1.5 GeV/fm 3 . 
The fireball's parameters discussed above are obviously in an interesting region. 
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2.2. N c and dN c /d?? Distributions 

Pseudo-rapidity distributions of charged particles have been measured in Au+Au [|J 
and Pb+Pb || collisions over nearly the full solid angle. The total charged particle 
multiplicity in central collisions increases from about 450 at AGS energy to about 1500 
at SPS energy. At AGS energy there are about equal numbers of pions and nucleons, 
while the pion to nucleon ratio is about 6 to 1 at SPS energy. This implies (within a 
thermal model, see below) a significant increase of the entropy per baryon in the fireball 
from about 15 to 38. 



3. Spectral Distributions 

Spectral distributions of high quality for protons and produced particles now exist |5|- 
[l~0"|j for heavy colliding systems both at AGS and SPS energy. In the following we will 
concentrate on the information on can glean from such spectra concerning stopping and, 
in particular, the presence or absence of collective features such as hydrodynamic flow. 

3.1. Rapidity Distributions and Baryon Stopping 

The rapidity distribution of identified baryons has been measured over the full solid 
angle at AGS energy ||[7|. The distribution is approximately Gaussian in shape with a 
peak at central rapidity and a width significantly narrower than that observed for the 
system Si+Al, implying strong baryon stopping. A similar difference is observed for the 
net proton rapidity distributions in Pb+Pb and S+S collisions at SPS energy Q. To 
get a complete picture here one has to await measurements of the A and A rapidity 
distributions. Current results from Na49 [|ll| are still somewhat controversial, with the 
A distribution nearly as narrow in rapidity as the A distribution. In any case, there is 
significantly increased stopping also at SPS energies, which is, of course, reflected in the 
increased entropy per baryon (see above). 

3.2. Transverse Momentum Distributions 

A general feature which has emerged from many measurements of transverse mo- 
mentum distributions is that, especially for central collisions of heavy nuclei, the in- 
variant distributions d 2 N / mtdm t dy which are approximately exponential in shape, i.e. 
d 2 N/m t dm t dy oc exp(— m t /T), have inverse slope constants T increasing linearly with 
the mass of the particle under consideration. Since the transverse momentum p t = m^ff3 t 
this fact has been widely interpreted as evidence for collective transverse flow: if there 
is a common flow velocity (or velocity profile) (3t superimposed on the random (thermal) 
motion of particles, the slope constant T, which is proportional to (pt) increases linearly 
with m. One can, however, also reproduce many features of the measured spectra by as- 
suming that the incoming nucleons undergo initial state scattering [Q. Since the amount 
of initial state scattering only depends on the number of nucleon-nucleon collisions each 
nucleon undergoes during the collision, hence is geometrical in nature, one does not need 
to invoke any collective expansion to explain the data in this approach. To resolve this 
issue, we discuss in some detail the two models and their implications. 
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3.2.1. Random Walk 

The aim of the random walk model |12| is to provide a description of minimum bias 
data on transverse momentum distributions in proton-nucleus collisions by a superposi- 
tion of the "kicks" the projectile nucleon undergoes in each nucleon-nucleon collision in 
the target nucleus. Since the extrapolation to AA collisions is then entirely geometric, 
one can test the predictions of this model directly by comparison to measured transverse 
momentum distributions. For pA collisions the approach is as follows. Each projectile 
nucleon undergoes iV^ nucleon-nucleon collisions in the target. Since, averaged over im- 
pact parameter 6, the nuclear thickness of the target is A/3R A , where R A = 1.12A 1 / 3 fm 
is the radius of the target, 

N A = 4/3R A a NN n . (1) 

Here n = 0.17/fm 3 is the nuclear density and the total nucleon-nucleon cross section ctjvjv 
is about 40 mb. One further assumes that each NN collision produces a fireball and that all 
final particles are emitted from the sequence of fireballs so created. Each fireball moves 
at a certain rapidity y and transverse momentum p t , determined by a boost invariant 
longitudinal expansion scenario (for y) and by the history of NN "kicks" for p t . To obtain 
the observed scaling with mass of the inverse slope constants the random walk proponents 
assume that the "kick" happens in transverse rapidity p t = 0.51n((m t + pt)/(rn t — Pt)) 
rather than in p t . Although this physically seems not obvious, it clearly introduces a 
"transverse velocity" distribution and, consequently, a mass dependent slope parameter. 
Since the sequence of kicks is stochastic, the distribution of the fireballs in transverse 
rapidity is assumed to be 

f pA {p t ) (xexp{-p 2 t /5l A ), (2) 
with 

5 2 pA = (N A - 1)5 2 . (3) 

Here, So is the only free parameter, apart from the fireball temperature (see below). It 
is determined by a fit to pA data and then kept constant for AB collisions. 

Generalizing the above expressions to nucleus-nucleus collisions we replace S pA by 

*AB(b) = (N c (b)/A + N c (b)/B - 2)51 (4) 
with 

N c (b)=T AB (b)a NN . (5) 

Here, the total number of collisions N c (b) for a given impact parameter b is evaluated 
with the help of the nuclear thickness function |U| . The purely geometrical nature of the 
problem becomes therefore transparent. Apart from a normalization factor which does 
not concern us here the expression for transverse mass spectra then reads: 

dN AB f 

mtdmt K j d Ptew(-Pt/$AB)g( m t), (6) 
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Figure 1. 

Comparison of (+) - (-) spectra for Pb+Au collisions with predictions of the random walk 
model for different centralities (increasing centrality from bottom to top). The data are 
from the CERES collaboration jlj . 
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Here, Jo and K\ are the modified Bessel functions and the particles are detected at rapidity 

V ru- 
in the following, we will compare predictions according to Eq. (||) with recent data on 
"proton" transverse momentum distributions, i.e. the difference of the spectral distribu- 
tions of positively and negatively charged particles, from the CERES collaboration [|14] 



These data for Pb+Au collisions were recorded with a multiplicity trigger corresponding 
to the upper 35 % of the geometrical cross section. This multiplicity range was subdivided 
into 7 exclusive intervals corresponding to (in a geometrical interpretation) mean impact 
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parameters of 2, 2.5, 3.5, 4.3, 5.1, 6.0, and 7.4 fm, respectively. The data are shown in Fig. 
|l].Data in different multiplicity (or impact parameter) bins were multiplied by consecutive 
factors of 8. The lines in Fig. [I] correspond to predictions of the random walk model ac- 
cording to Eq. (|6]), with 5q = 0.146 and a fireball temperature T = 150 MeV. Since we are 
only interested in the shape of the distributions here, the normalization of the calculation 
relative to the data is arbitrary. Note that, for central collisions corresponding to mean 
impact parameters 3.5 fm and less, the calculated distributions reproduce the shape of 
the data very well in the range 0.7 < m t — m < 2 GeV. For more peripheral collisions, 
however, the predictions by the random walk model exhibit inverse slope constants which 
are considerably smaller than seen in the data. It is important to realize that there is no 
parameter to change here, as the centrality dependence is determined by the thickness of 
matter traversed. If one determines the inverse slope constants from a fit to both the data 
and the calculations in the range 0.6 < m t — m < 2 GeV one can make this comparison 
more quantitative, as shown in Fig. Q. Here, the discrepancy between data and calcula- 




Figure 2. 

Comparison of measured inverse slope constants and values calculated in the random walk 
model for different centralities. Data are from the CERES collaboration fl4| . 
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tions becomes very obvious: certainly in its present form the random walk prescription is 
not consistent with the observed very weak impact parameter dependence of inverse slope 
constants measured for protons in Pb+Au collisions at SPS energy. Taking into account 
further that the calculated spectra also deviate strongly from the data at lower m t — m 
values we conclude that the random walk model is not anymore a contender to describe 
the flow-like features observed in the transverse momentum spectra. 

3.2.2. Transverse Flow 

For central collisions, the transverse momentum spectra of all observed particles can 
be well described in a hydrodynamical approach. The basic equations for azimuthally 
symmetric flow, i.e. central collisions, are very similar to Eq. |6|, except for the distri- 
butions in transverse rapidity. In a hydrodynamical approach, these distributions are 
obtained from the transverse flow velocity profile which, of course, depends on the initial 
conditions as well as on the underlying equation of state. We illustrate the success of this 
approach for data from the NA49 collaboration |]TTJ in Figs. ^] and [|. The hydrodynamic 



calculations used the transverse velocity profile recently calculated by Alam et al. [|T^ 
The initial conditions were chosen to reflect the measured charged particle multiplicities 
in central Pb+Pb collisions at SPS energy. For the equation of state it was assumed that 
a pure quark-gluon plasma undergoes a first order phase transition at T c = 160 MeV to 
a hadronic resonance gas which freezes out at constant temperature Tf. In Fig. |3] the 
data are compared to the calculations represented by the dashed lines for Tf = 120 MeV. 
Excellent agreement between the shape of the measured distributions and the calculations 
is obtained for all particle species ranging from pions to deuterons. We also show, in Fig. 
|3], the result of calculations with a linear profile tanhp t = (3 t = P ma x-^- Again, good 
agreement between data and calculations is obtained for Tf = 120 MeV and f3 max = 0.6. 
As is well known Tf cannot easily be determined from such fits, as higher Tf values can be 
traded off against lower (3 t values. This is illustrated in Fig. ^ where the solid lines show 
the results of a calculation with Tf = 140 MeV and (3 t = 0.45. We note, however, that the 
full hydrodynamic calculation does not exhibit this degeneracy: for Tf = 140 MeV, there 
is, with the present initial conditions and equation of state, no possibility to describe all 
measured distributions simultaneously. Under the assumptions discussed above, the data 
seem to clearly favor a relatively low freeze-out temperature near 120 MeV. The reason is 
that, to build up the flow in the hadronic phase, one needs the freeze-out to happen not 
too close to T c . A similar conclusion was recently obtained by the NA49 collaboration 
from a simultaneous analysis of tranverse momentum spectra and two-pion interferometry 



data ||16|| . 



3.3. Azimuthal Anisotropy and Directed Flow 

The study of collective flow effects via the experimental observation of anisotropics in 
azimuthal distributions for non-central collisions is motivated by theoretical predictions 
that a quark-gluon phase will lead to distinct differences in the flow patterns compared to 
what is expected for a hadron gas. Such phenomena were originally studied experimentally 



at Bevalac energies [17] and below. With the advent of the heaviest beams at the AGS 



and SPS such effects were also established there [|18| , |T9|| , at first in distributions of global 
observables such at Et and later also in distributions of identified particles. It is now 
customary to extract information on flow by determination of the Fourier coefficients Vi 
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Figure 3. 

Comparison of hydrodynamic calculations assuming a freeze-out temperature of 120 MeV 
with transverse momentum spectra for various particles. The data are from the NA49 



collaboration [1 IJ and are arbitrarily renormalized for clarity of presentation. The solid 
lines are calculations assuming a linear velocity profile. The dashed lines use the velocity 
profile determined by the initial conditions and the equation of state. For details see text. 



of the azimuthal distributions F((p) integrated over a certain rapidity interval 

n 

F(0) = F o (l + ^2^cos(#)). (8) 



i=l 



The dipole coefficient v i is also denoted as 'directed flow' and is associated to the mean 
transverse momentum in the reaction plane (p x ) by 

vi = (Px)/(pt), (9) 

the quadrupole coefficient t>2 is also called 'elliptic flow'. The odd Fourier coefficients 
change sign at midrapidity, while the even coefficients are symmetric. A quantity that 
was used to characterize the directed flow is the change of (p x ) with rapidity in the vicinity 
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Figure 4. 

Same as Figure 3 but for a freeze-out temperature of 140 MeV. 



of mid-rapidity and indeed the quantity 

F y = d(p x )/dy 



(10) 



is scale invariant and should be independent of beam energy if there is no change in the 
physics. Fig. |^ displays the information on directed flow of protons in terms of this scale 



invariant variable as a function of the beam kinetic energy per nucleon [17]] . One can see 
a plateau in the beam energy range of a few hundred MeV to 1 GeV and a steep fall-off 
to much lower values at the AGS and again much lower values at the SPS. It should be 
noted that the baryon density reached in the collisions is maximal for the AGS regime 
while it is lower and of comparable magnitude for the high end of the Bevalac/SIS regime 
and the SPS. So, at a qualitative level, one can clearly state that there is a difference in 
the equation of state over the energy regime considered with an increasing softening at 
the higher energies. 

It was shown at this conference [2(| that the directed flow of light nuclei is larger 
than the proton flow and increases monotonically with particle mass. The pion and kaon 
flow on the other hand appear to be a complex superposition of Coulomb effects, strong 
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Figure 5. 

Energy dependence from SIS to SPS energies of the proton sidewards flow as a function 
of the beam kinetic energy per nucleon. 



final state interaction, and collective flow which is apparent in a distinctly different p t 



dependence of V\ (see |pi] , P0|1 ). 

At AGS energy it was shown that in a hadronic cascade code (RQMD) one can in- 
deed reproduce the overall magnitude of the proton directed flow if a repulsive mean field 
between the nucleons is introduced |TH|. The same calculation however fails completely 
to reproduce the shape of the p t dependence of the proton directed flow which is on the 
other hand well matched by a hydrodynamically inspired picture of a sideways moving 
and expanding fireball [|T8| , pTj| . 

The elliptic flow of nucleons, quantified by the coefficient t>2, shows an interesting depen- 
dence on beam energy. At lower beam energies it reflects the energy dependence of the 
nucleon-nucleon interaction which is attractive at very low energies and becomes repul- 
sive in the 50 to 100 MeV/nucleon range. In terms of elliptic flow this is reflected in a 
change of sign of t>2, i-e. a 90° change in orientation of the ellipse (positive sign: long axis 
in the reaction plane) and leads to the so-called squeeze-out in the Bevalac/SIS energy 
regime with dominant nucleon emission perpendicular to the reaction plane. At AGS 
energy the elliptic flow is small but now oriented along the reaction plane and at SPS 
energy the elliptic flow is significantly stronger but oriented the same way. Such behavior 
had been predicted by theory |22j as a consequence of the changing relative importance 
of shadowing (leading to out-of-plane squeeze-out) and of collective flow related to the 
pressure build-up early in the reaction. At AGS and higher energies this compression 
related elliptic flow appears to dominate. The data are by now quantitative enough that 
a theoretical analysis should be performed to deduce the early pressure build-up. 
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Figure 6. 

Energy dependence from SIS to SPS energies of the elliptic flow (second Fourier coeffi- 
cient). 



3.4. Thermal Equilibrium and Freeze-out 

At the time of the Quark Matter '96 conference one of us reviewed the evidence for 
chemical and thermal equilibrium among the hadrons produced in ultra-relativistic col- 
lisions between heavy nuclei ||23|| . The conclusion at that time was that many features 
of the data indeed imply that a large degree of chemical equilibration is reached both at 
AGS and SPS energy. Since that time more data have become available (essentially all at 
SPS energy) and we will provide a brief update on the question of chemical equilibrium. 
We would like to stress, however, that a final word on this important question can come 
only after the data have been consolidated and when data cover as much as possible of 
the full solid angle. 

In Table [I] we show an update of presently available particle ratios for Pb+Pb cen- 
tral collisions at SPS energy, along with predictions from a thermal model calculation 
P^D which is slightly refined compared to our previous model [|25|| . Fairly good agrement 
between data and calculations is obtained assuming full chemical (including strangeness) 
equilibration for temperatures of 160 and 175 MeV, corresponding baryon chemical poten- 
tials of 200 and 270 MeV, and no strangeness suppression. In fact, the more consolidated 
data available now favor the larger baryon chemical potential; we mainly show the cal- 
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culation with T = 160 MeV and //& = 200 MeV as a reference (which best described the 
preliminary data available at the time of Quark Matter '96). The systematic uncertainties 
in the particle ratios as witnessed by the differences, e.g., in the A/A or H~ /A ratios from 
NA49 and WA97 are presently still too large to make a final judgement. We, therefore, 
feel that the claim of evidence for partial strangeness equilibration from a recent analysis 
of NA49 data [25| is premature. 



Table 1 

Particle ratios calculated in 2 versions of a thermal model for temperatures of 160 and 
175 MeV, and baryon chemical potentials fib of 200 and 270 MeV, in comparison to 
experimental data (with statistical errors in parentheses) for central collisions of Pb+Pb 
at SPS energy. 



Particles 


Thermal Model 


Experimental Dai 


a 




1 


2 


exp. ratio 


Exp. 


y 


p-p / neg a 


0.17 


0.23 


0.23(3) 


NA49 


0.2-5.6 


7T~ /lT + 


1.04 


1.06 


1.10(5) 


NA49 


all 


P/V 


0.086 


0.050 


0.055(10) 


NA44 


2.3-2.9 


0.099 


0.076 


0.085(8) 


NA49 


2.5-3.3 


d/d 


6.7 -10- 3 


2.1 -lO" 3 


3.6(8) -10- 3 


NA44 


1.9-2.1 


K+/K" 


1.57 


1.94 


1.61(15) 
1.85(9) 


NA49 
NA44 


2.5-3.3 
2.4-3.5 


K°/tt- 


0.145 


0.138 


0.125(19) 


NA49 


all 


A/A a 
A/A b 


0.147 


0.115 


0.128(12) 


WA97 


2.3-3.4 


0.17 


0.15 


0.19(1) 


NA49 


2.6-3.8 


20/(7T+ + 7T - ) 


0.020 


0.020 


9.1(10) -10- 3 


NA49 


all 


S + /S _a 


0.255 


0.273 


0.266(28) 


WA97 


2.4-3.4 


~-/A a 


0.114 


0.101 


0.093(7) 




n 


~-/A b 


0.093 


0.084 


0.13(4) 


NA49 


2.0-2.6 


S+/A a 


0.198 


0.239 


0.195(23) 


WA97 


2.4-3.4 


Q + /Q- a 

(n+ + n-)/(E+ + E-) a 


0.46 
0.154 


0.68 
0.168 


0.46(15) 
0.195(28) 


» 

J) 


n 
n 



a No feeding from weak decays. 

b Feeding from weak decays included. 



Note that, as discussed above, there is now evidence that, at SPS energy, thermal 
freeze-out happens at smaller temperature than chemical freeze-out. 

In Fig. [7|we show, for SIS, AGS, and SPS energies, the freeze-out parameters deduced 
from the chemical analysis of particle yields. It is interesting to note that, at the higher 
energies, the freeze-out points approach the calculated phase boundary p5| while freeze- 
out parameters for experiments at SIS energies never come close to it. 
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Figure 7. 

Freeze-out parameters deduced from the hadro-chemical analysis of particle production 
yields at SIS, AGS, and SPS energies. At the higher energies the freeze-out points ap- 
proach the phase boundary. 



4. Leptonic Observables 

4.1. The dilepton continuum 

First results from a measurement of the dilepton continuum at low mass in Pb+Au 
collisions were shown at the Quark Matter '96 conference P7[. Meanwhile, the CERES 



collaboration has released the final data [281. Similar to what was observed for the S+Au 



data, the electron pair yield in the invariant mass range 0.2 — 2.0 GeV is enhanced by 
a factor of 3.5 ± 0A(stat.) ± 0.9(syst.) over what is expected from neutral meson decays 
extrapolating from nucleon-nucleon collisions. The enhancement increases strongly (ap- 
proximately quadratically) with charged particle multiplicity and, hence, centrality. Most 
notable is that the observed enhancement seems to be concentrated at low transverse mo- 
menta. This is shown in Fig. [8] where, for Pb+Au collisions, inclusive e + e~ pair transverse 
momentum spectra are presented for three different pair mass ranges. The spectra are 
normalized again to charged particle multiplicity. For pair masses less than 0.2 GeV the 
data agree, as expected, with the sum of the hadron decay contributions. For larger pair 
masses, especially visible in the mass range between 0.2 and 0.6 GeV, the enhancement 
is strongest at very low pair transverse momentum. 

These results corroborate earlier findings for S+Au collisions of a lepton pair yield which 
is significantly enhanced compared to expectations for neutral meson decays. The new 
results, i.e. the strong centrality dependence and the concentration of the enhancement 
at low pair transverse momenta, should help to distinguish between various models put 
forward for the interpretation of these data. 
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Figure 8. 

The dependence on pair transverse momentum of the inclusive lepton-pair continuum for 
different mass windows. The data are for Pb+Au collisions from the CERES collaboration 



4.2. Charmonium suppression 

The NA50 collaboration reported anomalous J/^> suppression in Pb+Pb collisions at 
SPS energy at the Quark Matter '96 conference |29j. In 1996 the collaboration increased 
their data sample on J/\l/ by about a factor of 5, and modified their apparatus to cover 
a much larger range in centrality. The preliminary results reported at this conference 
are shown in Fig. |9|. Shown here is the ratio of J/^> to Drell-Yan production versus the 
geometrical mean path length L of the J/\I/ or cc state traversing the target and projectile 
matter. The quantity L is, apart from very central collisions where L saturates, a fairly 
good measure of centrality but, more importantly, can be used to compare various systems 
ranging from p+A to Pb+Pb collisions. 

Prior to data on Pb+Pb collisions the ratio of J/ty to Drell-Yan was found to decrease 
exponentially with L, implying absorption of the J/\l/ or its precursory cc state in the 
nuclear material with an absorption cross section of 6.2 ±0.6 mb. The new data shown 
in Fig. |S] confirm the findings reported at Quark Matter '96 but now also demonstrate 
that, for peripheral Pb+Pb collisions i.e. those with L less than about 8 fm, the ratio 
of J/\l/ to Drell-Yan follows approximately (actually the slope seems somewhat steeper) 
the same exponential behavior also called 'normal nuclear absorption'. In fact, the data 
exhibit a rather strong additional suppression setting in at around L = 8 fm. How rapid 
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the onset of anomalous absorption is needs to be studied further. However, it is clear from 
the correlation between transverse energy and impact parameter or L that fluctuations in 
transverse energy lead to uncertainties of L (even for very small transverse energy bins) 
of the order of 1 fm and, hence, no "discontinuity" can be expected in the data over L 
ranges less than that. 



m 


100 


■+ 
1 


90 


°l 




cm" 


80 


SZ 


70 


D 




>- 

_|_ 


60 


Id 


50 


o 






40 


? 
\ 


30 






b 








1 
QQ 


20 



** 1 58 GeV/c ** MRS 43 

A p(200 GeV/c)-A (A = W,U) (NA38-RESCALED) 

■ 32 S(32 x 200 GeV/c) - U (NA38-RESCALED) 

• 208 Pb(208 x 1 58 GeV/c) - Pb (NA50) 




Bff(y)/ff(DY) = A exp(-p a obs L) 

A = 57 ± 4 

a obs = 6.2 ± 0.6 mb 

x 2 /ndf = 0.13 



10 

L (fm) 



Figure 9. 

Anomalous J/ty suppression as observed by the NA50 collaboration [[30]]. For more details 
see text. 



Nevertheless, the most exciting interpretation of the anomalous absorption is that, in 
the center of the hot and dense fireball formed in the collision, the initially hadronic matter 
is converted into bubbles of quark-gluon plasma where bound charmonium states cannot 
survive. Certainly the observations are consistent with this scenario (see, e.g., the talk 
by Kharzeev [[H] at this conference). Alternative explanations have focussed on a scheme 
in which hadronic comovers break up the J/^> (see, e.g., the talk by Vogt [[3^] at this 
conference). However, within the framework of such models it turns out rather difficult to 
explain the complete absorption curve. Furthermore, for the absorption picture to work 
the comover density must be unrealistically high (on the order of 1/fm 3 ). 
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5. Summary and Outlook 

Over the past year the field has seen significant advances. There is now strong evidence 
from many experiments that the initially very hot and dense fireball formed in ultra- 
relativistic nuclear collisions expands with a common flow velocity prior to freeze-out. 
This is based on the analysis of single particle spectra as well as of two particle correlations 
which we could not cover here for reasons of space. Other collective features like directed 
flow and elliptic flow have now been also demonstrated at SPS energy. Both at AGS and 
SPS energy thermal freeze-out occurs at a temperature close to 120 MeV. Chemical freeze- 
out seems to occur at AGS energy at nearly the same temperature, while at SPS energy 
the corresponding temperature is about 170 MeV: apparently chemical freeze-out occurs 
close to where one expects the phase boundary based on the predictions of lattice QCD 
and simple bag models. New results from the CERES collaboration have consolidated 
the picture on an enhanced low mass electron pair continuum, although a real distinction 
between the various explanations put forward for the nature of the enhancement will 
probably have to await new data with much improved resolution and statistics. The J/^f 
suppression observed by the Na50 collaboration has not found any convincing explanation 
in terms of conventional scenarios. Clearly an interesting "picture" of ultra-relativistic 
nucleus-nucleus collisions is beginning to emerge. Based on this we look into the future 
with great enthusiasm. 
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